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Abstract 



We propose a self-consistent Ornstein-Zernike approximation for studying 
the Edwards- Anderson spin glass model. By performing two Legendre trans- 
forms in replica space, we introduce a Gibbs free energy depending on both 
the magnetizations and the overlap order parameters. The correlation func- 
tions and the thermodynamics are then obtained from the solution of a set of 
coupled partial differential equations. The approximation becomes exact in 
the limit of infinite dimension and it provides a potential route for studying 
the stability of the high-temperature phase against replica-symmetry break- 
ing fluctuations in finite dimensions. As a first step, we present the numerical 
predictions for the freezing temperature and the zero-field thermodynamic 
properties above freezing as a function of dimensionality. 
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I. INTRODUCTION 



This paper is the third one in a series devoted to the application of the self-consistent 
Ornstein-Zernike approximation (SCOZA) to classical spin systems with quenched disorder. 
The SCOZA has been formulated some time ago as a theory for simple fluid and lattice- 
gas systems [p]]. It is based on an Ornstein-Zernike (OZ) approximation for the direct 
correlation function c(r) that, by construction, enforces consistency between the different 
routes that give the thermodynamic potentials in terms of the pair correlation functions. In 
the language of magnetic systems, this means that the same Gibbs free energy Q (m, f3) is 
obtained whether one integrates the susceptibility x with respect to the magnetization m or 
the internal energy U with respect to the inverse temperature (3 = 1/ksT. In the case of the 
spin 1/2 Ising model with ferromagnetic nearest-neighbor interactions, this self-consistency 
requirement is embodied in a diffusive- like partial differential equation with m and (3 as 
independent variables. This equation could be solved only recently for the various three- 
dimensional lattices 0, showing that the SCOZA predicts all thermodynamic and structural 
properties with great accuracy, even in the close vicinity of the critical point ||. As the 
spatial dimension d goes to infinity, the SCOZA becomes exact and identifies with mean- 
field theory. This approach can thus be viewed as an approximate but non-perturbative 
way of taking into account thermal fluctuations in finite dimensions. Although this is an 
OZ theory which does not handle correctly long-range critical fluctuations, the asymptotic 
critical exponents are non-classical. For the 3-d Ising model, they have the spherical-model 
values in a very narrow region above the critical temperature, but are much accurate along 
the magnetization curve {(3 = 0.35) H|J. 

In the preceding papers of this series we used the replica method to generalize this 
approach to disordered spin systems. This allowed us to obtain an accurate description of the 
dependence of the critical temperature on dilution in the 3-d site-diluted Ising model and 
to study the influence of the random-field distribution on the phase diagram of the random- 
field Ising model (RFIM) for d > 4 pi. In both cases, possible replica-symmetry breaking 
(RSB) effects were ignored. In this work, we present a first application of the SCOZA to 
the Edwards-Anderson (EA) Ising spin glass model [|7| whose low-temperature properties 
remain a subject of controversies after nearly twenty- five years of intense activity. Whereas 
Parisi's mean-field theory || is generally accepted as the exact solution of the infinite- 
ranged Sherrington-Kirkpatrick (SK) model PlJlOll, there is yet no consensus on whether 
or not the RSB scenario associated with the appearance of multiple equilibrium states at 
low temperature survives in realistic, short-ranged, finite-dimensional models (for recent 



reviews of this problem, see |TT|,[T2|). A relevant open question is the existence of a phase 



transition in presence of an external magnetic field. Whereas the scaling approach based 
on the droplet model |T3|-[l5| predicts that a non-zero magnetic field destroys the spin glass 
phase, the mean-field picture suggests that there is still an Almeida-Thouless (AT) line JTH 
separating the spin glass phase from the paramagnetic phase. 

Although the SCOZA is based on a simple OZ ansatz for the correlation functions (so that 
the critical exponent i] is zero, a rather crude approximation for the 3-d EA model according 
to numerical simulations |17]]), it may be sufficiently accurate to give useful indications about 
the stability of the high temperature phase against RSB fluctuations. It is worth stressing 
that the SCOZA becomes exact when d — > oo, so that one recovers in this limit the behavior 
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of the SK model. Our purpose is thus to use this non-perturbative approach to generalize 
the AT analysis to finite dimensions, i.e. to determine the locus of singularities of the spin 
glass susceptibility xsg working from the replica-symmetric region. This amounts to study 
the eigenvalues of the wave-vector-dependent inverse susceptibility matrix at k = 0. As 
the inverse susceptibility matrix is just the matrix of direct correlation functions generated 
by a Legendre transformed Gibbs free energy, our program is the following: derive the 
appropriate OZ equations, assume an OZ form for the direct correlation functions (which 
means, in the present case, truncate them at nearest-neighbor separation), invert the OZ 
equations, and use all existing self-consistency requirements and exact relations (for instance, 
the values of the pair distribution functions at zero separation) to derive a set of partial 
differential equations whose solution will provide the phase boundaries in zero field and 
the AT lines. (Note that the SCOZA can also be combined with the Parisi RSB scheme; 
the main obstacle is the tractability of the calculation that requires inverting ultrametric 
matrices as done in Ref. [0.) This program, however, is not fully completed in the present 
paper which is only devoted to the calculation of the freezing temperature Tf and of the 
thermodynamic properties at and above Tf as a function of dimensionality. The stability 
analysis, which requires a more difficult numerical computation, will be presented later. 
The paper is organized as follows. In section 2, we introduce the model and we perform 
the double Legendre transform that allows one to introduce the direct correlation functions. 
The OZ equations in replica space are then solved in the case of replica symmetry. In section 
3, we define an OZ approximation for the direct correlation functions which, at the lowest 
level, coincides with standard mean field theory. We also derive the exact core conditions 
and self-consistency relations. In section 4, this formalism is used to derive the SCOZA 
equations in the case of zero field. Numerical results are presented in section 5. 



II. THE REPLICA-SYMMETRIC ORNSTEIN-ZERNIKE EQUATIONS 

We consider an EA spin-glass Hamiltonian for N spins <7j = ±1 on a d-dimensional 
hypercubic lattice 

Ti. = - ^ JijOiOj - H^cTj (1) 

<ij> i 

where the first sum runs over all nearest neighbor pairs and the couplings are independent 
Gaussian variables with mean J and variance J 2 ; H is an external magnetic field. 

Introducing as usual n replicas af (a = 1, ...,n) of the original spin variables and per- 
forming the average over disorder, one finds that the quenched free energy is given by 

T=-- lim i[exp(-/3^„) - 1] = -\ lim -[Tr exp(-(3H n ) - 1] (2) 

where 7i n is a temperature-dependent effective Hamiltonian 

^ = -D J °E « + E <^ ( 3 ) 

<ij> a a,b i a 
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and the trace in Eq. (2) is taken over the erf's. It is also useful to introduce a fictitious field 
A which couples to the quantity — 1/2 £V J2 a b < 7 i< 7 i- When A > 0, k B TA can be interpreted 
as the variance of a Gaussian random field but, in the present work, we are mainly interested 



in the limit A = (see Ref. JL8] for a study of the SK model in the presence of a Gaussian 
random field) We treat if and A as sources that we extend to site- and replica-dependent 
values iff and A® b in order to generate the correlation functions. We thus consider the more 
general Hamiltonian 

n n = - E t J o E « + p j2 E °» a A -EE H >i + E ( 4 ) 

<ij> a a<b i a a<b 

where we have restricted Af 6 to a < b with the condition Af b = A\ a and omitted the constant 
contribution — nN[c/3J 2 /4 + A/2] of the diagonal terms (c is the coordination number of the 
lattice). The averaged magnetization m" and overlap parameters gf 6 are given by 

< = =< < >n (5) 



and 

ab 9T n a h , 

% = -^ZF =< a i a i >n (6) 
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respectively, where < ... > n denotes the replica thermal average. At the end, we shall take 
the limit iff — > if and A" 6 — > A. Then, in the limit n — > 0, one has mf — > m, and, if replica- 
symmetry holds, qf — > q (which is the also the usual Edwards-Anderson order parameter 
q E A when A = 0). Accordingly, one has 9(JF '/ TV) / dH = —m and d{T /N)/dA = \(q — 1) 
(the —1/2 comes from the constant contribution omitted in Eq. (4)). 

The second partial derivatives with respect to iff and Af 1 define the matrix of (con- 
nected) correlation functions in replica space 

c) 2 f 

13 = ~ dH a dH b 1 3 >n ~ <a * >n< °i >n (7) 

1 J 



s~iab,cd 



dAfdAf 

1 J 



--< ViViVjVj 



>n ~ < 



~- a ~-b 

CT,: 



^ , „cd 

> n < a j a j 



<b,c< d) 



and similarly for G^' c {T n = (3T n -, H = (3H, and A = (3 A). 

We now perform a double Legendre transform that takes the fields iff and Af fc into mf 
and qf , respectively. This defines the Gibbs free energy 

e n = ?n + E E H >i + E E A *^ ab ( 10 ) 

j a i a<b 
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which satisfies Hf = dQ n /dm^ and A^ b = dQ n /dq^ b . Accordingly, in the limit n — > 0, 
mf — > m, and q^ b — > q, the Gibbs free energy Q(m, q, (3) satisfies 

and 

d{G/N) 1 



dq 2 



A. (12) 



Q n is the generating functional of the direct correlation functions (or proper vertices in 
field-theoretical language) 

QCtb _ _ (13) 

13 dmfdm b 



C%' cd = n 9 . g n n . (a<b,c<d) (14) 
* 3 dqfdq cd 



and similarly for C^ ,c (Q n = f3Q n ). 

The matrix C is the inverse of the matrix G and this defines a set of "Ornstein-Zernike 
equations" . In the limit of uniform fields, one has in Fourier space 

C ac (k)G cb (k) + C , °' cd (k)G' cd ' 6 (k) = 5 atb (16) 



c<d 



and 



C ab ' e {k)G e ' cd (k) + C ab > ef (k)G ef ' cd (k) = 5 a , c 5 b)d . (17) 

e e<f 

These equations are the starting point of our study and the first task is to solve them, 
i.e. to express the propagators G's in terms of the direct correlations functions C"s (or, 
conversely, the C"s in terms of the G's) for general values of n and then take the limit 
n — > 0. This is easy when replica-symmetry holds since G and C have the same structure as 
the Hessian or stability matrix of the SK model which has been analyzed by de Almeida and 
Thouless [rjfl (on the other hand, the problem becomes highly non-trivial when continuous 



RSB a la Parisi occurs, since one has to cope with the inversion of ultrametric matrices , see. 
e.g. |I2]|). There are seven different types of matrix elements, G aa , G ab , G ab ' ab , G ab ' ac , G ab ' cd , 



G a,ab = G ab,a^ Qc,ab = Qab,c ^ resp Qaa ? c ab ,...etc) and the two matrices of order n(n + l)/2 
can be easily block-diagonalized. There are three subspaces of dimensions 2, 2(n — 1) and 
n(n — 3)/2 which, following the standard terminology |l2[ , |i~9|| , we call respectively L (for 
longitudinal), A (for anomalous) and R (for replicon). The block-diagonalized matrices then 
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contain a single 2x2 L-block, n — 1 identical 2x2 A-blocks, and n(n — 3)/2 identical R- 
eigenvalues. In the new representation, the Ornstein-Zernike equations, hereafter called the 
replica-symmetric Ornstein-Zernike (RSOZ) equations, break up into three sets of equations. 
The first set reads 



G L (k) = C L (k)- 1 (18) 



with 



_ C aa + (n - l)C ab (n - l)C a > ab + (n ~ 1} 2 (n - 2) (>' fec | 

= " I 2C a ' ab +{n- 2)C a > bc C ab > ab + 2{n - 2)C ab ' ac + (n ~ 2) 2 (n ~ 3) C ,ab - cc ' ' 

Similar expressions hold for the elements of G . The second set reads 



G A (k) = C A (k)- x (20) 



with 



,a _( C aa - C ab (n - l)(C a ' ab - C a ' bc ) . , , 

— I n—2 ^Qa,ab (Ja,bc^ fjab,ab _|_ ^ ^(Jab,ac — (ji 3)(7 a6 ' crf ' ~ ' 

~ A 

and similar expressions for G . The last equation corresponds to the replicon sector, 

G R (k) = C R (k)- 1 (22) 

with 

(jR £jab,ab c^^ab,ac _|_ ^jab,cd (23) 

and a similar expression for G R . 

In the limit n — > 0, Eq (19) readily yields 

= C aa - C ab (24a) 

= C a ' bc - C a ' ab (24b) 

C£ = 2{C a > ab - C a ' bc ) = -2C[ 2 (24c) 

(7 2 L 2 = c ab ' ab - AC ab ' ac + 3C ab ' cd , (24d) 

and similarly for the four elements of the matrix G . They can be deduced from the elements 

of C by using Eq. (18) (for simplicity, in these expressions and in the following, we keep 
the same notations for G L , C L , G aa , etc... although the limit n — > has been taken). 

When n — > 0, we see from Eqs. (19) and (21) that C — > C (resp. G — > G ). But 

~ A * L 

the difference C — C is of order n, which yields 

1 -L -A / C ab —l(J a > bc \ 

i™ n^= ~ = ^ = \ 2C a > bc - C a > ab C ab ^ 2 - \C ab < cd ) ^ 

t / r^ab 1 riahc \ T 

' ' y 2G a ' bc + G a ' ab Qab,ac _|_ 3^ab,cd J ' / 
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This set of equations allows to express all the correlation functions C"s in terms of the 
G's (and conversely) in the limit n — > 0. The values of the direct correlation functions at 
k = in the longitudinal sector are directly related to the various susceptibilities. One has 

d 2 {Q/N) dH ~ L 

-^^=d^ = C ^ k = 0) (26a) 

dmdq dq 2dm ^ ^ 

d 2 {Q/N) ISA \* L 

^^ = -2^ = -2^ (k = 0) - (26C) 

In particular, when if = A = 0, C*n(k = 0) = /^x" 1 ; the inverse of the magnetic sus- 
ceptibility, whose divergence signals the occurence of the paramagnetic to ferromagnetic 
transition. More generally, in the approximate theory that is discussed below, the vanishing 
of the determinant of C L defines spinodal lines. On the other hand, the stability limit of 

the replica-symmetric solution (i.e. the AT lines) is signaled by the vanishing of C R (k = 0). 
In zero fields, one has 

C R {k = 0) = = xsh (27) 

where xsg — V^X^jt^ a i a j > t — < °i >t< &j >r) 2 ]av, is the spin-glass susceptibility 

It is also useful to note the consequences of the symmetry properties of the Hamiltonian 
7i n described by Eq. (4) for the correlation functions. Ti n is invariant under the global trans- 
formation {erf}, {H?} — > {—erf}, {—iff} for all spins in all replicas while leaving {A" b } un- 
changed. In consequence, averaged quantities involving an odd number of spins are zero when 
H = m = (even when A ^ 0). In particular, G a ' ab (r) = G a > bc (r) = C a ' ab (r) = C a ' bc (r) = 0. 
This implies that Gf 2 ( r ) = ^2i( r ) = Cvzi r ) = ^2i( r ) = an d the RSOZ equations involving 
the 2-replicas and 4-replicas correlation functions, Eqs. (16) and (17), decouple. TC n is also 
invariant in the global transformation {of}, {Hf}, {A" b } — > {— cr?}, {—H°-}, {—A" 6 } for all 
spins in a single replica a. This implies that G ab (r) = G ab ' ac (r) = G ab,cd (r) = C ab (r) = 
C ab > ac (r) = C ab,cd (r) = in the paramagnetic phase when H = A = 0. Finally, when 
Jo = 0, 7i n is also invariant in the local transformation {erf}, {Hf} — ► {— erf}, {—Hf} for 
a single spin i in all replicas while leaving {A" 6 } unchanged. This implies that m = and 
that G aa (r) and G ab (r) are local when H = (even when A ^ 0). More precisely, one has 
G aa (r) = 5 r:0 and G ab (r) = q5 r , . This implies that x = /^(k = 0) = (3(1 -q) in this case, 
as is well known |10| (the first equality follows from the fact that Cf 2 (r) = C^r) = 0). 



III. ORNSTEIN-ZERNIKE APPROXIMATION AND EXACT RELATIONS 

So far, all equations are exact provided that replica symmetry is unbroken. We now 
introduce an OZ approximation for the direct correlation functions. Since the interactions 
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in the effective Hamiltonian are restricted to nearest-neighbors, we assume that the seven 
distinct C's are non zero only at r = and r = e, where e is a vector from the origin to one 
of its nearest-neighbors. On the other hand, their dependence on m, q and (5 is not given a 
priori. We thus write that 

C*(r) = cg(m, q, P)5 r , + c?(m, q, (3)5 r , e (28) 

where x = aa, ab, aab, etc. ..The whole problem lies in the determination of the c 's and the 
ci's. 

When Jo = J = (or at infinite temperature with H and A fixed), the correlation 
functions are local (i.e. the c±s are zero) and can be calculated by averaging directly over 
the disorder distribution. Hereafter, we call this system the reference system. In terms of m 
and q, one finds that G% f (r) = (1 -m 2 )5 r , , G? ef (r) = (q-m 2 )5 r , , G"f/(r) = m(l-q)6 Tt0 , 

O r ) = (t-m)8 Tj0 , G^;\v) = (l-q 2 )5 r , , Gf e f c (r) = q(l-q)6 Tj0 , G?>f{r) = (r-g 2 )5 r , , 
with 



m 



/+oo 
Vx tanh(7f re/ + xK l J e )) (29a) 
-oo 



/+oo 
Vxttmh 2 (H ref + xK l J e )) (29b) 
-oo 

/+oo 
Vx tanh 3 (i7 re/ + xA^J) (29c) 
-oo 



/. 



+oo 

Vx tanh 4 (iJ re/ + xK l J e )) . (29d) 



■oo 



where = l/y^(27r)(ia;exp(— x 2 /2). Here, i? re / and A re / must be considered as auxiliary 
field variables that can be eliminated from Eqs. (29a) and (29b) to express t and r as 
functions of m and q, for instance as infinite double series, 

t(m, q) = 3mq - 2m{m 2 + 3q 2 ) + 6mq{2m 2 + 3q 2 ) - 6m(m 4 + 10m 2 g 2 + 13g 4 ) + ... (30a) 

r(m, q) = 3q 2 - 8q 3 - 2(m 4 - 16g 4 ) + 24g(m 4 - 7g 4 ) + ... (30b) 

The corresponding direct correlation functions are obtained from the RSOZ equations, Eqs. 
(20-27) (inverting the role of G and C). Their expressions are given in the Appendix. 

Now, the simplest approximation for the direct correlation functions at finite temperature 
is the Random Phase Approximation (RPA). It consists in setting the ci's equal to —f3 times 
the corresponding interactions in the Hamiltonian. Here, this gives 

C%>A<M) = cgTre/K 9, P) ~ ^A(k) (31) 
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Cg£(k) = c$? f (m, q, (3) - J 2 A(k) (32) 

where J = c/3J , J = c 1//2 /3J, and A(k) = l/c^ e exp(zk.e) is the characteristic function 
of the lattice. All other C"s remain equal to the corresponding C's of the reference system, 
CrpaO*-) = cg fe re/ (m,g,/3), C%p A (k) = c Q f ef {m,q, /?), etc... 
Using Eqs. (Al) and (A2), one then finds 

W k ) = (1 - 4g + 3r)D re/ - J A(k) (33a) 
C^Wk) = C 1 - - J 2 A(k) (33b) 

CS^paW = -\c^ RPA (k) =(m- t)D ref (33c) 

so 

dH 

— = (1 - Aq + 3r)D ref - J (34a) 
<9A 

^- = (1 - g),D re/ - J 2 (34b) 

W = -2^ = (m - t)0 "'' < 34C » 

It is straightforward to check that these results are also obtained by differentiating the 
following expression of the free energy 

TjN = (1 - qf + — m 2 - / Vx ln[2cosh(# + J m + x(A + J 2 q) 1/2 ] (35) 

where m and g are given by Eqs. (29a) and (29b) with H re f replaced by H + Join and A re / 
by A + J 2 q (the same modification must be performed, of course, in the definitions of t and 
r). Eq. (35) is just the replica-symmetric solution of the SK model in the presence of a 
Gaussian random field [T8|]. Therefore, the RPA is equivalent to mean-field theory when the 
Gibbs free energy is obtained by integration of the H- or A-susceptibility. (This statement 
is actually valid irrespective of the assumption of replica symmetry.) In particular, the 
instability with respect to RSB is given by the vanishing of 

Cgpx(k = 0) = C£,[l " 4 PA \ (36) 

where 

z« PA = (l-2q + r)J 2 . (37) 
When A = 0, the condition z p \ PA = 1 yields the usual AT lines | 16fl . 
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Mean field theory becomes exact when d — > oo. On the other hand, as is well-known 
in another context (see e.g., ||20||), the RPA in finite dimensions predicts pair distribution 
functions that do not satisfy the proper sum rules in the "core", i.e. when r = 0. Here, 
the exact core conditions are obtained by using the hard spin condition Oi = ±1 and the 
definition of q. One finds in the replica-symmetric case 

G aa {r = 0) = 1 - m 2 (38a) 

G ab (r = 0) = q-m 2 (38b) 

G a ' ab {r = 0) = m(l - q) (38c) 

G ab ' ab {r = 0) = l-q 2 (38d) 

G ab ' ac {r = 0) = q(l - q) . (38e) 

It can be checked that these relations are not satisfied by the RPA. This disease may be 
cured by adjusting the values of the corresponding direct correlation functions at r = 0, 
i.e. the values of Cq" 1 , Cq 6 , Cq'" 6 , c^' ab and c^' ac . In liquid-state theory, such an approximation 
is called the Optimized Random Phase Approximation (ORPA) J21J and is closely related 
to the mean-spherical approximation. We have seen in Ref. in the case of the RFIM 
that going from the RPA to the ORPA improves the predictions for non-universal quantities 
such as the critical temperature. It also modifies the critical behavior because of a subtle 



interplay between small-r and large-r correlations [22]. On the other hand, the ORPA 
(like the RPA) is not thermodynamically self-consistent as different Gibbs free energies are 
obtained depending on the route that is chosen for calculating Q from the pair correlation 
functions. In order to get a self-consistent theory we have also to adjust the values of the 
direct correlation functions at r = e, i.e. the values of the ci's. 

In the EA model, self-consistency is embodied in three Maxwell relations that can be 
obtained by considering the variation of Q as one varies the control variables J, m and q 
independently while keeping the ratio Jo/ J fixed. From Eqs. (4) and (10), one finds in the 
replica-symmetric n — > limit 

^§1 = -\&G»{r = e) + m 2 ) - (G^(r = e) + q 2 - 1)] . (39) 
dJ 2 4 J 2 

Then, using Eqs. (26) for the second partial derivatives of Q with respect to m and q, 
one finds that the cross-derivatives satisfy 



0Cft(k = O) J 1 d 2 r J, 



^G aa (r = e) - G abAb (Y = e)] (40a) 



dJ 2 2 J 2 Adm 2L J2 



0Cfe(k = O) 1 d 2 r J 



^G aa (r = e) -G ab ' ab (r = e)] (40b) 



dJ 2 Admdq'ji 
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^(k-0) = _i + I|l[i G -( r = e ) - G ab ' ab (r = e)] . (40c) 

Therefore, one has only five core conditions and three self-consistency relations available, 
whereas, from the OZ approximation, Eq. (28), there are forteen unknown functions c 's 
and Ci's to determine. In order to solve the problem completely, one must thus introduce 
additional approximations. 



IV. THE CASE H = J = 

As a first requirement, a sensible theory for the spin glass transition should yield rea- 
sonable predictions for the freezing temperature Tj in zero field. We thus take H = in 
the following, and, to keep things simple, we only consider the case Jo = 0. Then m = 
and the RSOZ equations, Eqs. (16) and (17), decouple, as noted earlier. Moreover, only the 
correlation functions depending on four replica indices are relevant, so that only two core 
conditions, Eqs. (38d) and (38e), and one self-consistency relation, Eq. (40c), come into 
play. In order to have only three unknown state-dependent functions to determine, we shall 
assume that cf ' ac (q, 0) = cf ' cd (q, 0) = and cf cd (q, 0) = cff e d f (q). According to the ORPA 
philosophy, these two additional approximations seem quite natural since only two distinct 
replicas interact in the Hamiltonian 7i n and there is no core condition associated to G ab,cd . 
Our simplified OZ approximation for the direct correlation functions in Fourier space thus 
reads 

C ab ' ab (k) = c^ ab (q,0) + 4 b ' ab (q,P)\(k) 
C ab ' ac (k) = c ab ' ac {q,(3) 

Since the direct correlation functions have the same spatial structure as in the RPA, this 
theory will reduce to mean field theory and become exact when d — > oo. 

It is now easy to calculate G ab > ab (r) and G ab ' ac (r) so to use the two core conditions and 
derive the SCOZA partial differential equation (PDE). Introducing the auxiliary function 
G D (r) = G ab ' ac (r) - 3/2 G ab > cd (r), one first notes that 

G aM6 (r) = -2[G£ 2 (r) - -G R (r) + G D (r)} (42a) 



G #c (r) = -\[G L 22 {v) - G R (v) + ^G D (r)} (42b) 
From Eqs. (41), one has 

C£(k) = C L 22 [l-^A(k)] (43) 

i riL ab,ab A ab,ac . ab,cd n ab,ab mL o- i i . • /^S^ 1 r*^ 

where Cg^ 22 = c — 4c + oc re j and z = —c x /G^ 22 . bmce the matrices Cj and C 
are diagonal when m — 0, one gets immediately 
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G&(r) = G L 22 P(t,z) (44) 
where G$ 22 = 1/Cq 22 and P(r, z) is the lattice Green's function defined by [23[ 



i /-7T ik.r 

P(r, z) = -. — r-T / d d k — — s . (45) 

1 ' (27r) d J_ n 1 _ z A(k) 1 ' 

Similarly, one has 

C R (k) = C*[l - * B A(k)] (46) 
with C,f = c[f ,afe — 2cQ b ' ac + and zr = — c" ' ob /C R . Hence, from Eq. (22), 

G R (r) = G*P(r,z R ) (47) 
with = l/Qf = (zr/ z)Gq 22 . Finally, from Eq. (25), one obtains 

with Gff = [3 cgg/2 - c o6 ' oc ](G£ 22 ) 2 = [egg + (1 - z/^)/G£ 22 ](G£ 22 ) 2 /2. Hence, 

G D (r)=G° j-izP(r,z)} (49) 

The core conditions, Eqs. (38d) and (38e), can then be used to express Gq 22 and Gq in 
terms of z and z R . One finds 

L _ 2z(l - q) 
G °< 22 " 3z R P(z R ) - zP(z) m 



and 



G ° " *P'(*) + P(z) [1 ~ g - 2Z 3z R P(z R )-zP(z) ] (51) 
where P(z) = P(r = 0, 2). Moreover, z and z B are related via 

(1 - q)[3z R P(z R ) - zP(z)] 2 = 2 z[3^P(z fl ) - «P(*)M*) + (1 - -)^-{zP{z))\ 

z R dz 

+ ^l-q)^(zP(z))cf;f f ( q ) (52) 

The self-consistency relation, Eq. (40c), can now be expressed as a PDE in the unknown 
function z(q, A), where A = J 2 = c(3 2 J 2 . One obtains 

^(k = 0) _ l5 V^r-e) f53l 
dA - _1 "2^ G (r " e) (53) 

with 
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AWl. n\ 3z R P(z R ) - zP(z) 

^ 22 (k = 0) = ^— ^ (1 - *) (54) 



and 



Gab , ab(r _ e) _ (1 _ _ U2 agM - gg(g) - 1 m 2 , g(f) i g(f) 1 

(55) 

where we have used that P(r = e, 2) = [P(z) — l]/z. 

Given the appropriate boundary conditions, the set of Eqs. (52-55) allows to calculate 
z and zr and thus the correlation functions and the Gibbs free energy for all values of q 
and A. However, from the definition of the Green's function, Eq. (45), one must have 
< z < 1 and < Zr < 1. In particular, z = zr = corresponds to the high-field (A — > 00, 
q — > 1) or high-temperature (A — > 0) limit, while z = 1 and zr = 1 define respectively the 
spinodal line Qf 2 (k = 0) = and the limit of stability of the replica-symmetric solution 
in the q — A plane. At the freezing temperature Tf when A = 0, one has simultaneously 
z(q = 0) = zn(q = 0) = 1. Indeed, c^^(g = 0) = 0, so that 2: = zr is always solution of Eq. 
(52) when q = (this is not the only solution, but the only one that is physically acceptable). 
Accordingly, G R (q = 0) = Gg 22 (g = 0) so that c^ ,cd (q = 0) = 0, as required by symmetry 
(cf. the end of section 2). Another exact requirement is that A(q = 0) = for every A in the 
high-temperature phase. Since A re f(q = 0) = and dA/dX = — q — (l/2)dG ab ' ab (r = e)/dq 
(the integrated form of Eq. (53)), this implies that 



8 

£tab,ab i 

dq 



= (56) 

<?=o 



in the high-temperature phase. We take this equation as the boundary condition for the PDE, 
Eq. (53), on the line q = 0. 

This suggests to take G ab ' ab (r = e) as the unknown function in the numerical integration 
of Eq. (53). The PDE is then rewritten as a non-linear diffusion equation, 

dG f|jj 1 d 2 G 

d\=-im {1 + 2W ) (57) 

oz I q 1 

where C = (k = 0) and G = Q ab > ab ( r = e). The integration is carried out by a simple 
explicit algorithm in which the "space-like" variable q and the "time-like variable" A are 
discretized, and the partial derivatives are approximated by finite difference representations 
The first derivative with respect to A is used to update G at the temperature step n+1 
by evaluating the second derivative with respect to q at the step n. For given values of q 
and G, z and zr are obtained from Eqs. (52) and (55) by using a standard Newton-Raphson 
algorithm (note that these equations does not depend explicitly on A and that c^ e Jq) can 

be tabulated once and for all). C^fe = 0) is then obtained from Eqs. (54). Although 
an implicit method, such as the one used in Ref. |3[], would have the advantage of being 
unconditionally stable, we have achieved numerical stability above and at Tf by using small 
grid spacings, Ag = 5.10 -3 and AA = 10 -5 . 
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V. RESULTS AND DISCUSSION 



As one decreases the temperature (i.e., increases A) in the numerical integration of the 
PDE, there is a value of A for which z(q = 0) = Zn(q = 0) = 1. This corresponds to the 
divergence of the spin glass susceptibility, and it defines the freezing temperature Tf. The 
SCOZA results for Tf are compared in Table 1 with the predictions of a 1/d series expansion 



|25| for d = 9, 7, 5 and with a recent Monte Carlo estimate [26] for d = 3 (we have only 
considered odd dimensions because the Green's function P(z) is easier to evaluate). One 
notes a significant improvement over the predictions of mean- field theory (T / m/ = y/2d). In 
fact, the results are comparable to those of the Bethe-Peierls approximation |27| for d > 5 
and better for d — 3. Note again that the present theory is exact in the replica- symmetric 
high-temperature phase when d — > oo and thus predicts the correct Tf in this limit. It is 
likely that the small deviations for d = 9 and 7 come from the fact that only c1 b,ab (q,/3) 
is nonzero in the simplified OZ approximation, Eq. (41). Indeed, the exact calculation of 
(jab,ab^ r ^(jab,ac^ anc j (jab,cd^ ^ order 1/d shows that the three functions do extend to 
nearest-neighbor separation PS|. But the introduction of c" ,ac and c^ ,cd would require self- 



consistency relations involving the three direct correlation functions separately. This is only 
possible if one breaks replica symmetry in some way. Below the upper critical dimension 
d — 6, the fact that rj is no longer zero comes also into play and may explain the more 
significant deviations from the simulation results, especially for d = 3 (rj ~ —0.5 according 
to Ref. |I7||). 

When J = 0, the zero-field internal energy per spin is given by f3u = —(A/2) [1 — 
G ab ' ab (r = e)]. Therefore, at the freezing temperature, Eq. (55) readily gives Uf/J = 
—cJ/(2P(l)kBTf). The corresponding numerical values are also given in Table 1. Finally, 
as an illustration of the behavior of the specific heat above Tf, we plot c(T) = du/dT for 
d = 5 and 3 in Fig. 1. Note the broad maximum that appears above the freezing temperature 
for d = 3, in agreement with experimental and simulation data [[K|]. This is in contrast with 
the mean-field behavior that predicts a cusp at Tf. 

In summary, the SCOZA results for the properties of the EA model at and above the 
freezing temperature as a function of dimensionality are in reasonable agreement with the 
known estimates. This represents an encouraging step for proceeding further in studying 
the stability of the replica-symmetric solution at temperatures below Tf. However, we have 
not yet been able to obtain a reliable solution of the SCOZA PDE below freezing and this 
work is still in progress. 



14 



APPENDIX A: DIRECT CORRELATION FUNCTIONS IN THE REFERENCE 

SYSTEM 



From the expressions of the G's in the reference system and Eqs. (18) and (24) , one 
finds 

Cft,™/ = (1 " 4fl + 3r)£W (Ala) 

C22,ref = (1 - Q)Dref (Alb) 

= -\ckref = (m- t)D ref (Ale) 

where 

D re f = C'll,re/C , 22,re/ + ^(Pl2,ref) 2 

= [(1 -g)(l -4g + 3r) + 2(m-t) 2 ]~ 1 (A2) 
and r and £ are given in Eqs. (29). Moreover, Eq. (25) yields 

C tref ~\ C tref \ = n2 / 1 ~ 4g + 3r m-t 



r, a,fec a,af> afe,ac 3 ab,cd I ^re/ 1 _0( m _ f) \ — n 

ZC 0,ref ~ C 0,ref C ref ~ 2 C 0,re/ / V Z v" L ) 1 1 

m 2 — q ^(t — mq) \ / 1 — 4g + 3r m — t 

m(l + q)-2t -q + |(3r - g 2 ) J ' \ -2(m - t) 1-g 



(A3) 



Finally, Eq. (22) for the replicon sector gives 

65/ - T^Tr < A4 > 

It is now the straightforward to calculate all the C re /'s from these equations. The results 



are 



c o"e/ = DrefK 1 ~ 4g + 3r) 2 (l + m 2 - 2q) + 2(m - t)(l - 4g + 3r)(2m + mq - St) 

+ (m-£) 2 (6g + g 2 - 6r - 1)] (A5) 



= Ce/ " (1 " 4« + 3r)^«/ (A6) 



CoSf = 2{(1 - 4g + 3r)[(m(q - m 2 )(l - g) - i(i -mq)(l + q- 2m 2 )] 

2 



g 2 3r 



+ (1 - q)(m - t)(q + ^ - —) + (m - t) 2 (2t - m - mq)}D 2 ref (A7) 
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Colfef = C tef + (* " "0 A-e/ (A8) 

o 2 3r 

= 2{(1 - g) 2 (g + y - y ) + (m - t) 2 (l + 2m 2 - 3g) + 2(1 - q)(m - t)(2t - m - mq)}D 2 ref 

+ 2 (l-g)(^-g) + (^-^ D (A9) 
l-2g + r re/ V ; 



ab,ac _ ab,cd , (jj-jQ^-g) + (m - t)' 
C 0,re/ - C 0,re/ + ! _ 2 g + r re/ 1 UJ 



afe,ab _ ab,ac ab,cd 1 / A i i \ 

C 0,re/ - ZC 0,re/ C 0,ref + 1 _ 2 + _ V ' ' 
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d 


k B T s f coza /J 


r , 1 

k B T b f e /J 


uf oza /J 


9 


3.91 


3.82 a 


-2.16 


7 


3.35 


3.22 a 


-1.91 


5 


2.66 


2.41 a 


-1.63 


3 


1.58 


0.95±0.04 6 


-1.25 



Table 1: Freezing temperature and internal energy per spin at Tj of the d-dimensional 
EA spin glass model with Gaussian couplings and J = 0. The SCOZA results for Tj are 
compared with the best series or simulation estimates, a: Ref. [p5 |. b: Ref. |26]. For 
d = 9,7, 5, Tf is obtained from the 1/d series estimates of the freezing temperature for the 
±J distribution using the transformation formula given in Ref. [25 |. 
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FIG. 1. SCOZA predictions for the specific heat vs reduced temperature T/Tj of the EA spin 
dass model with Gaussian couplings and Jq = for T > Tf. 
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